Interplay between positive feedbacks in the generalized CEV process 
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The dynamics of the generalized CEV process dXt = aX^dt + bX^dWt (gCEV) is due to an 
interplay of two feedback mechanisms: State-to-Drift and State-to-Diffusion, whose degrees are n 
and m respectively. We particularly show that the gCEV, in which both feedback mechanisms are 
POSITIVE, i.e. n,m > 1, admits a stationary probability distribution P provided that n < 2m — 1, 
which asymptotically decays as a power law P{x) ^ with tail exponent ji — 2m > 2. Furthermore 
the power spectral density obeys S{f) ^ jj, where /3 = 2 — 2(l^-i) ^ ^ > 0- Bursting behavior of the 
gCEV is investigated numerically. Burst intensity S and burst duration T are shown to be related 
by 5'-T^ 



The dynamics of the state Xt of a system which 
is open to a rapidly fluctuating environment can 
be described by the non-linear stochastic differen- 
tial equation 

dXt = f{Xt)dt + g{Xt)dWu (1) 

Wt the standard Wiener process, under the assump- 
tion that 1.) noise enters linear, and 2.) the White 
Noise approximation is valid, see The drift 

and the diffusion 'coefficients' are depending on the 
recent state Xt and hence represent 'State-to-Drift' 
or 'State-to-Diffusion' feedbacks, respectively. The 
resulting dynamics and consequently properties 
like the stationary pdf of the gCEV, the spectral 
density, and also burst statistics are shown to be 
due to the interplay between these two feedback 
mechanisms. 

The following (informal) argument shows that if 
both feedback mechanisms f{X) and g{X) have a 



particular functional relation to each other, given 

by 

f{X) = ag{X)g'{X) (2) 

then their interaction generates a power-law like sta- 
tionary probability distribution - if it exists - in that 



(Here and in the following the notion F{x) ~ f{x) 
means that the function F{x) = f{x) for large x.) 
Note that the proportionality factor a enters the co- 
efficient of the power-law tail. The process consid- 
ered in dXt = b'^ (m - |) X^'^'^dt + bX^dWt, 
corresponds to g{X) = and a = 1 — ^ so that 
the stationary pdf decays as a power-law according 
to P(x)-^. 

I. THE GENERALIZED CEV PROCESS 
(GCEV) 
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In the following we consider a particular setting 
which is that the drift and diffusion coefficients obey 
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f{Xt) = a and g{Xt) = bX^^. In this case one 
obtains the Ito diffusion process with positive drift 
and diffusion parameters a, b given by 

dXt = aX^dt + bX^dWt, n,m>0. (4) 

This process is a generahzation of the Constant 
- Elasticity - of - Variance model, dXt = 
aXtdt^bX]^ dWt^ which was originahy proposed by 
Cox, Ingersoh and Ross to describe the dynamics of 
interest rates in an equihbrium economy and which 
plays an important role in Mathematical Finance, 
see references below. The gCEV process, eqn|4] de- 
scribes dynamics by the superposition of two differ- 
ent feedback scenarios: One is the State-to-Drift 
feedback incorporated in the deterministic part of 
the gCEV 

(5) 

while the other one concerns State-to-Diffusion 
feedback due to 

dXt = bXJ^dWt (6) 

In the following, we will focus on the case that both 
dynamical components exhibits positive feedback 
simultaneously in that we require 

1 < n,m, < oo! (7) 

In this case gCEV dynamics results from the in- 
terplay of two POSITIVE feedback scenarios, each of 
which generates self- amplifying, i.e. 'explosive' be- 
havior in itself. This is easily seen that the drift 
term with positive feedback gives rise to a Finite- 
Time-Singularity, i.e. ^ oo as t ^ tc, where 
tc = 7^-^o~"^7 while the solution on the finite in- 
terval [0,tc] is Xt oc 1 . Positive feedback in 

the state-to-diffusion term also leads to bursting be- 
havior, in that Xt can attain arbitrary large values 
while it always remains finite. This follows from the 



fact that the solution of eqn [6] is the inverse power 
of a d-dimensional Bessel process Xt oc ^ — , 

iiBiir"' 

where B is a d dimensional Brownian motion, where 
d = > 2, for details see [5 . Since therefore 

d{m) > 2, it follows from the transitivity of B^, that 
B^ escapes to oo for t ^ oo slower than t a.s., see 
[1 , while the origin is polar, i.e. it will not be 
touched by B^. On the other hand B^ has a positive 
probability to visit any finite e-neighborhood of the 
origin before escaping. Consequently the dynamics 
exhibits arbitrary high but finite excursions. 
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FIG. 1. Time trail of the CEV process dX ^hX^dW 

Hence, for n, m > 1, both singularities are en- 
tirely different: While positive feedback in the state- 
to-drift component leads to a 'real' Finite-Time- 
Singularity in that ^ oo within [0, tc\^ Xt remains 
finite even if feedback in state-to-diffusion is posi- 
tive. Nonetheless, if both positive feedback play in 
concert, the process exhibits a fat-tailed stationary 
probability distribution, provided that the State-to- 
Diffusion feedback is positive (m > 1) and strong 
enough with respect to the State-to-Drift feedback, 
i.e. m > ^(n + 1). Bursting behavior is reflected 
in that the stationary pdf decays as a power law 
P{x) ~ ^ for large x with an exponent obeying 
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jii = 2m > 2. For large the tail exponent only 
depends on the state-to- diffusion feedback parame- 
ter m, while the state-to -drift feedback parameter n 
determines the pdf only for small x. 

II. STATIONARY PDF OF A 
GENERALIZED CEV PROCESS 

A. The CEV process 

The standard CEV process is obtained from eqn [4] 
for n = 1 

dXt = aXtdt + bX^dWt, a,6>0, (8) 

A typical time series generated by the CEV process 
for m = I is shown in Figjl] An extensive discus- 
sion of the CEV process and its relation to other 
processes including Bessel processes can be found in 
[U |3] . While most of former research on the CEV 
process has been restricted to the case < m < 1, 
we instead focus on the case that state-to-diffusion 
is subject to positive feedback m > 1. A detailed 
discussion of this process and the following theo- 
rem as well as the proof can be found in [5]. As 
shown there, a CEV process with m > 1 is equiva- 
lent to a radial Ornstein-Uhlenbeck process for order 
u(m) — o/ IN and hence admits a closed form an- 

V f 2(m — 1) 

alytical solution given below: 

Result 1 (Solution of the CEV process for 

m > 1). The unique and strong solution of the CIR- 
CE V model, eq^ with m > 1 is 

1 



X{t) = c(m) 



|l/(m-l)' 



(9) 



\mt) 

with c{m) = ^6(m — 1)^ ^ , where M(t) he 
a d- dimensional mean-reverting Ornstein-Uhlenbeck 
process, whose dimension is a function of the feed- 
hack parameter m given by 

1 



Its components obey d Mi (t) = —jiMidt^dBiit) with 
a > and B{t) the standard Wiener process, while 
its square norm is ||M(t)||2. 

The proof is based on the observation that the 
Lamberti transform of this process [S] takes the form 
of a radial Ornstein Uhlenbeck process of order 
v{m) = 2{ra-i) ' i^^- -'-^ ^^^^ noie we show that the 
CEV process eqn. [8]with positive state-to- diffusion 
feedback (m > 1) admits a stationary probability 
distribution, which is uni-modal and asymptotically 
decays as a power-law with its tail exponent propor- 
tional to m only. 

Result 2 (Stationary pdf for the CEV process 

for m > 1). Let dXt = aXtdt + bX^dWt [CEV] be 
defined on the non-negative reals [0, oo) with a, 6 > 
and Wt the Standard Wiener process. Then, if m > 
1 a stationary probability distribution exists and is 
similar (not equal) to a Type-2 Gumbel distribution 



P{x) 



where c 



2a 



-2(m-l) 



2(m-l) 



(11) 



b^2i^) >0 andM= ^^^^ <oo IS 
a normalization constant with 7 = > 1 • The 

stationary pdf takes its unique maximum in x^ = 
( - — ^ ) where ly(rri) = 1^ is the index of 

\ a m—l J V ^ 2(m— 1) 

radial Ornstein Uhlenbeck process equivalent to eqn. 

m 



Proof. As shown in [5^ the solution of the CEV pro- 
cess [CEV] is an inverse power of a radial Ornstein 



Uhlenbeck process of dimension d = 

for 1 < m < 00, given by Xt cx ^ —j— 



> 2 
whose 
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S{m) = 2 



m — l 



> 2 



(10) 



components obey dMi = —aMidt + dWi. Since 
< ||M||2 < 00, and 00 are natural boundaries 
for the CEV process with m > 1, thus the probabil- 
ity current over these boundaries is zero. □ 

For m = I the stationary probability distribution 
asymptotically decays as a power law P{x) ~ -\ 
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FIG. 2. The stationary pdf belonging to the Ito process 
dX = aXdt + bX'^dW for m = 3/2 (red) and m = 2 



(green), see eqn. 11 defined on the non-negative reals. 



for large more precisely 



P{x) 
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The pdf is shown in Fig. [3j 



(3 1,2 X 



(12) 



P(x) 



10"^ 
10"' 



10 



10 



-12 . 



16 




10"^ 10^ 10^ 10^ 10^ lO"^ 10^ 10^ 



FIG. 3. The stationary pdf belonging to the CEV pro- 
cess dX — aXdt + bxidW, with positive feedback 
on 'state-to-diffusion. Simulated data in green, while 
the red dashed curve is due to eqn. |12| 

The equivalence between the CEV model with 



m > 1 and the rOU process of index ^{m) = 2(^-1) 
shows also up in the functional form of the station- 
ary pdf. In fact the distribution can be rewritten 
as 



P{x) OC X 



-(2+l/z/(m)) 



exp 



-u{m) X 



which shows that the exponential part of the pdf is 
controlled by the order of the corresponding rOU 
process. Furthermore, the stationary probability 
distribution belonging to the CEV process with m > 
1 is uni-modal for all b and decays for large X as a 
power-law, i.e. P{x) ~ Its graph is sketched in 
Fig [2] As an immediate consequence we have 

Result 3. Given the CEV process with feedback pa- 
rameter m > 1 and let v{m) = ttt^:^— rr- Then 

1 

where e > is a small constant. Consequently 
—00 < P{m) < 2 for 1 < m < 00. 

This follows from the more general case, see below. 
Note that for m = 3/2 and c = 0, /3 = 1, so that in 
this case S{f) ~ 4 



S{f) 



I3{m) 



2(m-l) • 



(13) 
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B. The generalized CEV process with positive 
feedback n, m > 1 

The interplay between state-to-price feedback and 
state-to-diffusion becomes obvious when consider- 
ing the Fokker-Planck equation belonging to the 
gCEV process. Note that by transforming Xf into 



Yt 



6(m-l) X 

a{Yt)^dWt with drift 



T, the gCEV process becomes dYf 



a(^b{m-l)Yt^ 



(14) 



m-lYt_ 

The corresponding potential U{y) = — a{y')dy' 
of the corresponding Fokker-Planck equation is of 
the form 

1 . 
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FIG. 4. The power spectral density of the CEV process 
dX — aXdt + bxidW. Simulated density in red, while 
the green dashed lines is the 1// decay. 

where rj = > if n < 2m — 1 and positive 

otherwise. Thus there is a bifurcation occurring at 

= 2m — 1, such that U{y) is convex for n < 
and concave for n > n^, see Fig[5j Particularly, for 
n < n*, and oo are repelling and the potential has 
a unique minimum. On the other hand, for n > n*, 
both and oo are attracting, while U{y) ^ ^ 
— oo with some positive c for y approaching 0. In 
terms of X this means that for n > n^^ that small 

^ 0, while large ^ oo. 

While for n < the gCEV process can be re- 
garded as a diffusion trapped in a convex potential, 
i.e. behaves locally similar to an Ornstein-Uhlenbeck 
process, one can expect that in this case it admits a 
stationary probability density. 

Result 4. The generalized CEV process admits a 
stationary probability distribution P{x) ifm > 1 and 
n < 2m — 1. The stationary probability distribution 
yields P{x) = M exp (-c where 
c = I? 2m-n-i ^' '^^^^ normalization constant 
^ l^^ifey^ ^ (2m-i)-n ^^^^ asymptoti- 




FIG. 5. Potential U{ii) of the Lamperti transformed 
gCEV process; the upper curve corresponds to n < n*, 
while the lower one is for n > n* . 

cally decays as a power law with tail exponent /i > 2 



Note that this excludes the Geometric Brownian 
Motion case [n = m = 1], while it includes the 
case that the process exhibits positive feedback 
of both: state-to-drift feedback [n > 1] as well as 
state-to-diffusion [m > 1]. In fact, given the degree 
n > 1 of positive state-to-drift feedback, then the 
degree of state-to-diffusion feedback m must be 
positive as well while sufficiently large m > ^(n — 1). 

Fig [g] shows the empirical distribution P (red) 
compared to the analytical solution, see eqn [TT] 
(green dashed lines) for the case m = 3/2 and 
n = 4/3, i.e. for the gCEV process in which both 
partial processes, eqnn [5] and |6] exhibit positive feed- 
back. 

The power spectrum density of a gCEV asymp- 
totically decays as a power-law with tail index /3, 
which is a direct function of the index u of the re- 
lated Bessel process and is, particularly independent 
of the feedback exponent n of the drift term. 
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Note that for m = | and c = 0, the power- 
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FIG. 6. The stationary pdf belonging to the gCEV pro- 

4 3 

cess dX — aX^dt + bX'^dW, with positive feedback 
on both 'state-to-drift' and 'state-to-diffusion. 



spectrum shows pure 1/f behavior, see Fig. [7| 
while for m = |, the spectruni is flat, f3 = 0. The 
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Result 5. The gCVE process eqn with feedback 
parameters m > l,n < 2m — 1 admits a power spec- 
trum 

S{f) ^ j^, 13 = 2 - {l + e)v{m), (15) 

where e > and v{m) = 2{m-i) ^"^ index of 
the radial Ornstein-Uhlenbeck process equivalent to 
the CEV process with m > 1 and c > is a small 
parameter. 

Proof. The proof follows from a results in 16] by not- 
ing that the gCEV process for large Xt can be approx- 
imated by dXt = c{X)X^'^-^dt + bXJJ'dWt, where 
forn < 2m+l the coefficient c{X) = :^X^~^'^'^^^-^ 
is almost constant for large Xt and approaches 
from above if Xt oo. Therefore approximating 
the g CEV process for large but finite X by 



dXt 



cX^'^-^dt 



bXJ^dWt, oO (16) 



we obtain eq. (3) in 16] with the substitution c = 
(m — I). According to eqn. (33) in ^6] the spectral 
density reads S{f) ^ , /3 = 1 -\- 2(^-1) • Inserting 
v = 2{m — c) we obtain that for the process in eqn 



16 (3{m) = 2 



l+2c 
2(m-l) 



which gives the result eqn 



putting v{m) = 2(^-1) ^ — 



FIG. 7. The power spectral density S(f) ^ ji belonging 
to the gCEV process dX = aX^dt + bX^dW. 

n-dependence of the spectral density shows up only 
for small /. One can show that S{f) for small / is 
an increasing function of n. 



III. BURSTS GENERATED IN THE GCEV 
PROCESS 

Regarding the transformed gCEV process Yt for 
n < as a diffusion being trapped in a convex 
potential U (y) as in Fig. [s] makes clear that the 
dynamics of Xt allows for a sequence of arbitrary 
high but finite outbursts even on short time scales, in 
agreement with Fig[l] Since the bursting behavior, 
i.e. Xt large, is governed by the state-to-diffusion 
feedback parameter m, we can restrict ourselves to 
the case n = 1 for investigating statistical properties 
of burst. That is, we will numerically consider the 
CEV process 

dXt = aXtdt^bx} dW 
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in the following. Kaulakys and Alaburda [4 consid- 
ered the case m = 2 in dX = p- x^^~^dt + x^dW, 
for X > Xm > 0. Note that is the case X is 
distributed according to a power law with tail 
exponent 2 (m — p-), as follows from eqnjsj In this 
particular setting they found numerical evidences 
for clear power-law statistics of bursts. Since in our 
case, power-law behavior only exists asymptotically, 
we can expect power-law burst statistics only 
asymptotically. 



FIG. 8. Intermediary power-law regime in PxiT) for 
different threshold values: x = 2 (red circles), x = 4 
(green squares) , x = 8 (blue triangles) . 

The size of a burst is defined as 5 = dt, i.e. 

the integral over the super-threshold trajectory in 
the burst period T. It turns out to be related to the 
burst duration T by 

5 - (17) 



A BURST is regarded as a super-threshold event: 
Let {Xt) be the solution of a gCEV process. The 
hurst interval of the k-th burst is defined as the 
time interval between crossing the threshold x > 
from below and the smallest time at which the 
threshold is crossed back from above. By a slight 
abuse of notation we also denote the length of this 
burst period by T/^. In Fig |8] its probability distribu- 
tions P^(T) = P{Tk = T\x) are shown for different 
threshold values: red x = 2, green x = 4, blue x = 8. 
The distribution P{T) of burst durations admits an 
intermediary power-law regime with -\. 
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FIG. 9. Size 5* of a burst and its duration T for three 
thresholds, x = 2 (red circles), x = 4 (green squares), 
X = 8 (blue triangles) . 
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